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In their recent Letter Wang and Landau have pro- 
posed a very effective MC method of producing evenly 
distributed histograms for the degeneracy 5 (£'i) of the en- 
ergy levels of different discrete spin models. The method 
suffers, however, from a drawback. In the course of their 
calculations Wang and Landau have to reduce their mul- 
tiplication factor / to a value f final which differs from 
1 only by 10~^. At this stage the algorithm comes very 
close to entropic sampling. Indeed when g{Ei) has the 
value 10^ then the modified density of states (DOS) 
g{Ei) — *■ g{Ei) * f final in the new algorithm and the 
modified DOS: g{Ei) — > g{Ei) -f 1 in entropic sampling 
are identical. Thus the new method has the definite ad- 
vantage of rapidly building a rough estimator for the DOS 
but in the end the work of refining it to useful accuracy 
is as time consuming as entropic sampling. 

We have ourselves used a similar method as Ref. |^ 
to calculate histograms of the degeneracy of g{Ei) of 
three dimensional Ising models and of the degeneracy 
g{Eo,Mi) of 2d-Ising models for a fixed value Eq of 
the energy as a function of the magnetization Af^. We 
have also produced histograms for the three dimensional 
ANNNI-model. All these data are unpublished. 

In our approach the acceptance probability of a spin 
flip is also governed by the DOS which has been accu- 
mulated so far. In this respect our approach is identical 
to that of Wang and Landau, but it differs from Ref. 
in two important aspects: (i) the histogram g{Ei) of 
the DOS is updated (multiplied by /) only every tenth 
step, (ii) the transition observables T{Ei,A) are also 
evaluated by recording all attempts to move from an en- 
ergy Ei to Ei+ A with A = 0, ±4, ±8 for the 2d-Ising 
model. On taking the logarithm one obtains the entropy 
S{Ei) = In g{Ei) from the DOS and the inverse tem- 
perature 13 = ln [T{E, + A, -A)/T{Ei, A)] /A from the 
transition observables. It turns out that the latter data 
suffer much less from statistical errors than the former. 
In contrast to Ref. ^ we do not vary / - it is kept con- 
stant during the simulation. In the examples which follow 
we have added 77 = 0.004 viz. rj' — 0.002 to the entropy 
S{Ei) after every 10 attempted spin flips. This leads to 
multiplication factors f = oi 1.004 viz. 1.002. 

It turns out that the error one makes is roughly propor- 
tional to T], i.e., \ngrj{Ei) « S{Ei) + rif{Ei) + Kj^ where 
r]f{Ei) is this error and Kr/ is an unimportant constant. 
A very good approximation for S{Ei) can therefore be de- 
termined from two MC runs with different values of rj and 
7?': SiE,) « [r]lngr,'{Ei)-rj'\ngrjiE,)]/ii^-r,') + K. 


Fig. 1 has been produced for the 32 x 32 2d-Ising model. 
It shows the second derivative d^s{e)/de^ of our simu- 
lated data. Here s = S/N and e = E/N are the entropy 
and the energy per spin. These data scatter around the 
exact data of Beale [Q. The difference to these data is 
shown in the inset. Gaussian smoothing over 20 channels 
reduces the maximum deviation to 0.005. 
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FIG. 1. The second derivative d?s{e)/de'^ and the differ- 
ence to the exact data (inset). 

As the entropy is determined up to an unimportant 
constant, the derivatives of the entropy are more as- 
sertive. Therefore one should not compare the entropy 
itself, but its derivatives with the exact results |2|. We 
have chosen the second derivative for this comparison be- 
cause it enters into the microcanonically defined specific 
heat c{e) = —s''^{e)/s"{e) and it is thus very important 
that high quality data are obtained for it. 

In summary: Because of our substraction scheme we 
can keep the entropy increment and concomitantly the 
effectivity of the simulation procedure five orders of mag- 
nitude above Ref. 

We thank Michael Promberger for his advice. 
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